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We investigate the influence of a shear flow on the process of nucleation. Mesoscopic nonequilib- 
rium thermodynamics is used to derive the Fokker-Planck equation governing the evolution of the 
cluster size distribution in a metastable phase subjected to a stationary flow. The presence of the 
flow manifests itself in the expression for the effective diffusion coefficient of a cluster and introduces 
modifications in the nucleation rate. The impfications of these results in condensation and polymer 
crystallization are discussed. 

I. INTRODUCTION 

In Paper Tm of this series we used mesoscopic nonequilibrium thermodynamics (MNET) to analyze the effect of the 
inhomogeneities of the medium on nucleation. Another interesting case in which the medium may exert an influence 
on the nucleation process arises from the presence of flows or stresses in the system. This factor is specially relevant 
in crystallization, which often involves mechanical processing of the melt, such as extrusion, shearing or injection 
that may affect the crystallization process drastically. Therefore it is desirable to implement a model that takes these 
mechanical influences into account in a correct and consistent way. Several models have been proposed to capture 
some features of phase transitions under shear flowi but until now, there has been no theory providing a complete 
description of crystallization under these conditions. 

The objective in this paper will be to analyze precisely the possible effects of mechanical stress in nucleation and 
crystallization, focusing on the simplest influence. To that end, we discuss the case of nucleation in the presence of a 
stationary flow, which for simplicity is chosen to be a shear flow. MNET>2i4 will be used to derive the kinetic equation 
governing the dynamics of nucleation in the presence of a shear flow, whose influence on vapor-phase nucleation and 
polymer crystallization will be evaluated qualitatively. 

It is worth pointing out that the presence of the flow destroys the isotropy of the system, and leads to a distinction 
between diffusion in different directions. Moreover, it introduces spatial inhomogeneities, which induce spatial, and 
velocity and temperature fluxes. Coupling between these currents may be very important, and consequently the 
evolution of the probability density, velocity and temperature fields will be governed by a highly coupled set of kinetic 
differential equations. These equations can be derived within the general framework of MNET^^, following the steps 
developed in paper I of this series^. For illustrative purposes, it is best to restrict the analysis to effects originating 
purely from the presence of the flow. Consequently, for the sake of simplicity, we assume isothermal conditions. 

The paper is organized as follows. In Section II, we derive the Fokker-Planck equation governing the evolution 
of the size and velocity distributions of the emerging clusters in the presence of stationary shear flow. In Section 
III, we simplify the description, focusing on the diffusion regime. We thus obtain the kinetic equation governing 
nucleation in the presence of a stationary flow. The implications of our results for nucleation experiments and in 
polymer crystallization are discussed briefly in Section IV. Finally, we summarize our main conclusions in Section V. 

II. NUCLEATION IN A STATIONARY FLOW 

Consider a metastable phase in which the emerging clusters (liquid droplets or crystallites) are embedded. Assume 
that this phase acts as a heat bath imposing a constant temperature on the system. Following the scheme developed 
in Paper li, the goal of this section is the derivation of the kinetic equation that describes nucleation in the presence of 
a shear flow. For this purpose, we will work within the framework of MNET. To account for spatial inhomogeneities, 
we perform a local description in terms of /(n, x, u, t) /N, the probability density of finding a cluster of n G {n,n + dn) 
molecules at position x S (x, x -I- c?x), with velocity u G (u, u + du), at time t. 

The starting point is the Gibbs equation for the entropy variation of this system with respect to a reference state. 
In this situation, the reference state is a stationary one characterized by a steady shear fiow velocity profile vo(x). 

Assuming local equilibrium, the entropy variation is given by the Gibbs equation (see Appendix I) as 
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TpSs — pSe + ppSp ^^Pj l^^cdndxi. (1) 

Here s(x,i) and e(x,i) are the local entropy and energy per cluster, respectively, p{x,t) is the hydrostatic pressure, 
p{x,t) is the total number of clusters per unit of volume, x,u, t) is a non-equilibrium chemical potential, and 
c(n, X, u, t) — f{n,x,u,t)/p is the number fraction of clusters. Eq. (QJ remains vaHd within an element of mass 
following the center of mass motion of the cluster 'gas—. The explicit expression for the chemical potential can be 
obtained, as in paper I, by comparison of the Gibbs equation with the Gibbs' entropy postulate^, yielding 

f 

p{n,u,x,t) ^ ksTluj i-pieqix), (2) 

Jleq 

where pieq{x) denotes the local equihbrium chemical potential, which is independent of the n and u, and fieq is the 
local equilibrium distribution function corresponding to the reference state, a function given by 

r , ^ f pieq - C{n,u)\ 

Jleq (n, u, x) = exp ( I . (3) 

The quantity C(n, u) denotes the free energy of formation of a state described by the internal variables (n, u, ), and 
is given by 

C(n,u) = AG(n) + im(n)(u-vo)^ (4) 

where AG(n) represents again the energy of formation of a cluster of size n at rest and the second term is its kinetic 
energy with respect to the steady state velocity profile vo. 

The next step toward obtaining the Fokker-Planck equation, is the formulation of the conservation laws for the 'gas' 
of clusters. In the absence of external body forces, the distribution function obeys the following continuity equation 

which introduces the phase space current J^, arising from the interaction of the clusters with the surrounding 
metastable phase. 

In order to derive the entropy balance equation, we also need the expression for the variation of energy with respect 
to its value at the stationary state. The presence of the external fiow is responsible for the appearance in that equation 
of a term for 'viscous heating', that gives rise to variations in the temperature field. To maintain isothermal conditions, 
and following ideas introduced previously in the implementation of the so-called 'homogeneous sheari£*i, we assume 
the existence of a local heat source capable of removing heat generated in the process. Under this assumption, the 
energy of a volume element remains constant along its trajectory so that its balance equation can be omitted in the 
following analysis. 

Introducing into the Gibbs equation Q the equation for number fraction 

,_^_V./(u-vo)-^.J„-^, (6) 
derived in Appendix II, we obtain, after some straightforward algebraic manipulations, the entropy balance equation 

ds 

where the entropy flux Jg is given by 

Js^-ks J f {In f - l){u- \o)dndu- ^ J C{n,u)f{u-\ofdndu, (8) 
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and the entropy production, a, which must be positive according to the second law of thermodynamics, is 

a = [ 3u ■ ^dndu - ^ [ dnd\x 
1 J ou 1 J on 

(u - vo) • V QTO(n)(u - vo)2^ dndu. (9) 
In order to obtain the entropy flux we have additionally employed the identity 

J V • ((u - vo)/) In / dndu = V • ^(u - vo)/(ln / - 1) dndu-pS/ ■ vq , (10) 

and have omitted, in the entropy production, the term— (fcsp — y)V ■ vo arising from the bulk viscosity, which, in 
general, is negligibly small. 

The entropy production can be interpreted as a sum of conjugate forces-fluxes pairs. In the present situation, the 
forces are ~5^f^j and — ^^V (im(n)(u — vq)^) while their conjugated fluxes are the velocity flux J„, the flux 

in size space J„, and the relative spatial current 3x = /(u — vq), respectively. 

The next step is the formulation of coupled linear phenomenological equations relating these fluxes and forces. The 
situation under consideration presents some subtleties. Spatial and velocity currents have the same tensorial rank 
(they are both vectors) which implies that, in general, they must be coupled. Moreover, the presence of shear flow 
destroys the isotropy of the system by introducing a special direction. Consequently, in general, phenomenological 
coeflicients are expected to be tensors. In addition, locality in the internal space will again be assumed . Taking these 
considerations into account, the expressions for the currents are 

m „ \ 1 /. 

Jx = +j; jOxx -y^o ■ (u- Vo) + j; i--ux ■ (11) 

1 dfi ^ / X 

Jn = i--uu ■ + T^L-ux ■ Vvo • (u - Vq), (12) 

where we have used the Onsager reciprocal relation 

^UX — ^XU ; (■^^) 

as well as the result 

V Qm(u - vo)2^ = -mVvo • (u - vq). (15) 

It is useful to redeflne the phenomenological coefficients in a more convenient way. By identifying Z?„ = as 

the diffusion coefficient in rt-space, which in turn can be identified with fc+(n), the rate of attachment of monomers 

— ^ /■ — \ p 

to a cluster, while introducing the friction tensors a = and ( — "^p^" as well as the explicit form of the 

chemical potential, equation the currents J„ and J„ can be written as 



C - Vvo 



^dn ksT On 

Introduction of these expressions for the currents into the continuity equation yields 



u-vo -^^-Z, 16 

TO OU 
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C ■ Vvo 



dn 
/(u - vo 



-/ 



1 dC 

ksT dn 

^ksT df 
a TP 

TO OU 



(18) 



which is the Fokker-Planck equation describing the evolution of the inhomogeneous density distribution of clusters in 
the presence of steady flow. 



III. BALANCE EQUATIONS IN THE DIFFUSION REGIME 

The Fokker-Planck equation, Eq. Ijl8|l . we have derived retains information concerning the evolution of the cluster 
velocity distribution. Although the movement of clusters may play a significant role in the nucleation process, as 
discussed in Ref4, velocity dependences can hardly be measured because velocity distribution relaxes to equihbrium 
so rapidly. Therefore, as discussed in Paper I, a simplified description of the process can be provided, by coarse- 
graining the dependence on the velocity variables. The relevant quantities are then the different moments of the 
distribution function, namely the density of clusters 

fc{n,x,t) ^ J f{n,x,u,t)du , (19) 

the coarse-grained velocity of the clusters 

yr,{n,x,t) = f-' Jufdu, (20) 



and the second moment 



V = y"/(u-Ve)(u-Ve)ciu, (21) 

which is related to the kinetic definition of the pressure tensoi»Si2ii£. 

Proceeding along the lines of Section III of Paper I, integration of Eq. II18|I over velocities leads to the corresponding 
balance equations, namely 

^ = ■ fc^c- f Jndu , (22) 
dt dnj ' ^ ' 

for the density of clusters in n-space. 



fc^ - -V • T' - S • (Vc - vo) - y(u - Ve)^ J„du, 



(23) 



for the momentum, and 



/ (U - Vc)(u - Vc)^ J„rfu, 



(24) 



for the pressure tensor, where Q is, as in paper the kinetic part of the heat fiux and B = a + C ■ Vvq. 

A remarkable feature of these equations is the fact that the balance equation for the pressure tensor, Eq. 1I2H1 . 
is similar to the one obtained for the thirteen moment approximation to the Boltzmann equationpiS. It is thus 
worth pointing out that the mesoscopic nonequilibrium thermodynamic treatment of the problem has been able to 
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recover results similar to those of kinetic theory. In particular, inertial regimes and relaxation equations for the 
main quantities can be obtained within this framework. This result reinforces the fact that application of the well- 
established postulates of non-equilibrium thermodynamics as indicated in Refi^, suffices to provide a general scheme 
under which non-equilibrium processes of both macroscopic and mesoscopic nature can be treated. 

The second term on the right hand side of the momentum balance equation can be identified with the frictional 
force per unit mass exerted on the cluster suspended in the host fiuid. B plays the role of the friction tensor. The 
modulus of this tensor again determines the characteristic time scale for the relaxation of the velocities. In general, 
such a friction tensor can be calculated from hydrodynamics; however, for the sake of simplicity, we approximate this 

tensor by S ~ pt, the Stokes diagonal friction. This identification leads to a = /3l — C ■ Vvq. 

The inverse of the friction coefficient /?~^ again shows the separation of the dynamics of the system into a short-time 
inertial regime and the longer time diffusional regime. For times t > (3~^ the system enters the diffusional regime. In 
this situation, inertial effects manifested through the presence of the time derivative can be neglected. Moreover, 
the contribution arising from the current J„ in the velocity relaxation can be neglected since the rate of relaxation of 
velocities is usually faster than that of cluster sizes determined by D„. 

The results derived thus far are valid for any stationary ffow profile vq. However, to proceed further, and for the 
sake of simplicity, we will focus on the particular interesting case of a shear fiow. 

Proceeding along the lines indicated in Paper I and in Refpii, we find that, in the diffusion regime, the expression 
for the pressure tensor with a shear flow present, is 



(25) 



From this equation, we conclude that Brownian motion of the particles contributes to the total pressure tensor of the 
suspension in two ways. The flrst contribution is the well-known scalar kinetic pressure given by 



P = /c, 26 

m 

which is the equation of state for an ideal 'gas' of clusters. The second contribution comes from the irreversible part 
n° of the pressure tensor, a part that can be written in the form 



(1 + C ) • Vvo 



(27) 



where Dq = is the diffusion coefficient of the cluster for the liquid at rest, and the superscript indicates a 
symmetric traceless tensor. This last equation defines the Brownian viscosity tensor 



r^B = Dofc{t+ C), (28) 
which contains the "Brownian viscosity" Dof^, and the contribution due to coupling with the non-equilibrium bath, 
proportional to ( . 

Inserting the expression for the pressure tensor in the velocity balance equation ll23ll , and taking the diffusion limit 
for the cluster current we obtain 



(29) 



where 



V = Da 



(30) 



can be identified with the spatial diffusion coefficient, which in this situation possesses tensorial character. Using Eqs. 
H28II and Ij26|l . V can be rewritten as 



V = D. 



Vvo 



(31) 
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Inserting this expression into the equation for the density, Eq. H22|l . one finally obtains 

which is the kinetic equation for nucleation in presence of a shear flow in the diffusion regime, and where, as discussed 
in Paper I, AG{n) is the effective nucleation barrier in the diffusion regime given by 

AG(n) = AG(n) + ^fcsTlnn, (33) 



IV. INFLUENCE OF SHEAR FLOW IN NUCLEATION AND CRYSTALLIZATION 

In this section, we analyze the conditions under which mechanical stress that cause a shear flow in the system, 
may influence the nucleation process. Obviously, a shear rate modifles transport and results in alterations in density, 
pressure and temperature proflles. To isolate the influence of mechanical stress, heat and mass transfer will be 
neglected. The focus will be on the direct effect of flow on nucleation rate, in both condensation and crystaUization, 
with special emphasis on polymer crystallization. 

A. Influence on Condensation 

As shown in the previous section, an important influence of flow on nucleation is exercised through its effect on the 
diffusion coefficient as exhibited by Eq. 113 Typically, nucleation rates J are given hy^^^ 

AG("») 

J- /c+(n*)e"^?^, (34) 

where fc+(n) is the rate of attachment of molecules to a ti— cluster, and AG{n*) is the height of the nucleation barrier 
(i.e. the free energy of formation of the critical cluster of size n*). The significance of the modification of the diffusion 
coefficient lies in the fact that the rate of attachment of molecules to the nucleus, fc"*", and thus the pre-exponential 
factor in the nucleation rate, Eq. Ij34|l . is roughly proportional to the the diffusivity of moleculesi^. The fact that 
the diffusion coefficient could be altered by the fiow implies that the nucleation rate is modified accordingly. The 
modification, prescribed by Eq. l(3T|l . depends on the shear rate, the pressure and the viscosity, and its magnitude is 
expected to be important when 

-|VvoI>l, (35) 
P 

i.e., for high enough values of viscosity, and shear rate, and for low pressures. 

To be more precise, let us estimate the magnitude of the correction for a typical experiment in a laminar fiow diffusion 
cloud chambeuifi. Representative values of the shear rates, the pressure and the viscosity in these experiments are 
Vwo ~ 10'^ s~^, p ~ 10** dyn/cw?, and rj ~ IQ^^ poise. These numbers yield an approximate value of 10~^ for the 
order of magnitude of the correction introduced by the shear. Consequently, in this range of values, neither diffusion 
nor the nucleation rate are significantly altered by shear. 

B. Influence in Polymer Crystallization 

Whereas in condensation under normal conditions, the modification of the diffusion coefficient due to the shear 
fiow is in general neghgible, this is not the case in crystaUization and specially in polymer crystallization. The 
presence of fiow affects crystallization in different ways. On the one hand, nucleation rates can be directly modified 
by fiow. Typical values of shear and injection rates, and of the mechanical stress applied to a polymer melt during 
crystallization are usually rather large. In addition, the viscosity of the supercooled melt can also be very large. 
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specially near the glass transition. In this situation, the adjustment of the diffusion coefficient, given by Eq. (PH, can 
be very significant and thus modify the nucleation rate dramatically. 

On the other hand, the symmetry breaking inherent in the flow causes the diffusion coefficient to no longer be a 
scalar and to develop a tensorial nature. The rate of the process then depends on direction so that nucleation, and 
specially the crystal growth, are no longer isotropic. 

Finally, another important factor to be considered is that flow also induces stress on the surface of the cluster that 
may fragment droplets that become large enoughii*!^. This effect can be easily incorporated into our description by 
including in the free energy of cluster formation the additional energetic cost associated with the shear stress. 

V. CONCLUSIONS 

The main objective of this paper has been the analysis of the kinetics of nucleation and crystallization when the 
metastable phase is subjected to a flow resulting from application of mechanical stress. Using MNET, the kinetic 
equation that describes the evolution of the cluster size distribution in a stationary flow was derived. For simplicity, 
our results were speciallized to the case of a shear flow, and its influence on condensation and crystallization was 
evaluated. 

The main effects that a shear flow exerts on nucleation can be summarized as follows. First, such flow alters 
transport and consequently the evolution of the cluster distribution function, that among other things, determines 
nucleation and growth rates. Second, flow changes the cluster spatial diffusion coefficient in accordance with equation 
(jnOJ. Typical values of the parameters controlling this adjustment indicate that the effect is not very important in 
condensation. However, high viscosity and other peculiarities of polymer crystallization suggest that a shear may 
promote drastic changes in the process. The formalism developed in this paper can help to explain experimental 
results of polymer crystallization under sheaui^. However, quantitative prediction is limited, for the moment, by the 
unavailability of the magnitudes of the physical parameters required in calculation. Measurement of these parameters 
has to be performed before a complete characterization of the effect of flow, using our theory, on the kinetics of 
polymer crystallization can be explicitly evaluated. 

As indicated, the main object of this paper has been the analysis of the effects that flow by itself may induce on 
the nucleation process. For this purpose, isothermal conditions were assumed. A more complete and more realistic 
treatment that models mechanical processing of the melt may be carried out by incorporating thermal effects and 
their coupling to the shear flow, in the manner discussed in Paper I. 
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Appendix A 
Local Thermodynamic Relations 

The Gibbs equation for a multicomponent thermodynamic system is^ 



TSS = 5E+ pSV - fiSN. (36) 
In the continuous Umit this expression becomes 



TSS = 5E+ p5V - J ^iSN{n, u) dndu, 



(37) 



where different values of n and u denote the equivalent of different "species", and N{n, u)dndu represents the number 
of particles of a given "species". 

Now construct the local version of Eq. l(37ll . For this purpose, we will describe the system in terms of f{n, x, u, t), 
the number density of clusters of size n G {n,n + dn) at x G (x, x + rfx), having velocity u S (u, u + du), at time t. 
We define the specific thermodynamic quantities s, e, v, and c as follows 



/ 



spdx (38) 
epdx (39) 
V ^ / Idx = / vpdx (40) 



N{n, u) = y /o'x = J cpdx (41) 

where p(x, t) — J f dndu is the number density of clusters per unit volume at x, while v — p^^, and c — j;- 
Using these quantities, the local expression of the Gibbs equation becomes 

TS{ps) = 5{pe) + p5{l) - j p5f dndu, (42) 

which is the form used in Paper I. 

Alternatively, this equation can be written as 



(43) 



Tp6s — pSe — ppSp ^ ''^ P J P'^c.dndu + ^Ts ^ e — pp ^ + j iicdndu^ Sp = 0, 
which, using the definition of the Gibbs energy to cancel the last term, reduces to 

Tp6s = pSe + pp6p^^ ^ ^ J 1^^'^'^^^^' 
which is the expression employed in the present paper. 
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Appendix B 
Balance Equation for the Number Fraction 

Integration of Eq. © with respect to the cluster velocity u and size n, leads to the macroscopic equation of 
continuity, which can be written in the form 

| = -.V.v, (45, 

where p(x, t) is the density of the cluster 'gas', given by 

p(x, ^ J f dndu. (46) 
In Eq. H45ll.vfx. t) is the average velocity of the Brownian particles defined by the expression 

pv(x, i)= J uf dndu, (47) 

and where the "total derivative" is 

4 = -T- + V ■ V. (48) 

dt dt ^ ' 

Using the continuity equation, Eq. J^l, and Eq. II45II . one obtains the equation for the number fraction 

dc df f dp \ 9 1 ,Af^\ 

Under the assumption of local equilibrium , the average velocity v coincides with the velocity of the reference state, 
i.e. V = vq, thus leading to Eq. 
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